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Abstract 

In 1942 I. J. Schoenberg proved that a function is positive definite in 
the unit sphere if and only if this function is a positive linear combination 
of the Gegenbauer polynomials. In this paper we extend Schoenberg's 
theorem for multivariate Gegenbauer polynomials. This extension derives 
new positive semidefinite constraints for the distance distribution which 
can be applied for spherical codes. 



1 Introduction 

Let M be a metric space with a distance function r. A real continuous function 
fit) is said to be positive definite (p.d.) in M if for arbitrary points pi, . . . ,p r 
in M, real variables x\, . . . ,x r , and arbitrary r we have 

r r 

i=l j=l 

or equivalently, the matrix (/(%)) b 0, where the sign y stands for: "is 
positive semidefinite" . 

Let S" _1 denote the unit sphere in K™, and let ifij denote the angular dis- 
tance between points pi,Pj- Schoenberg [L7] proved that: 

/(cos if) is p.d. in S™ -1 if and only if f(t) = J2kLo fkG^\t) with all ff. > 0. 
Here G^\t) are the Gegenbauer polynomials. 

Schoenberg's theorem has been generalized by Bochner [5] to more general 
spaces. Namely, the following fact holds: / is p.d. in a 2-point-homogenous 
space M if and only if f(t) is a nonnegative linear combination of the zonal 
spherical functions (see details in |12., Th. 2], [51 Chapter 9]). 

Note that the Bochner - Schoenberg theorem is widely used in coding theory 
and discrete geometry for finding bounds for error-correcting codes, constant 
weight codes, spherical codes, sphere packings and other packing problems in 
2-point-homogeneous spaces (see [8j [12l Q31 HSJ [16] and many others) . 

Consider the following problem: Let Q = {qi, . . . , q m } be a set of points in 
M. To describe the class of continuous functions F(t,u,v) in 2m + 1 variables 
with t € M, u,v £ R m , F(t,u,v) — F(t,v,u) such that for arbitrary points 
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Px,...,p r in M the matrix (F(tij,Ui,Uj)) >: 0, where tij — r(pi,pj), U{ — 
(r(Pi,qi), ■ ■ .,T(pi,q m )). 

Denote this class by PD(M, Q). If Q = 0, then PD(M, Q) is the class of p.d. 
functions in M . 

Recently, Schrijver |18j improved some upper bounds on binary codes using 
semidefinite programming. Schrijver's method has been adapted for some non- 
binary codes (Gijswijt, Schrijver, and Tanaka [10]), and for spherical codes 
(Bachoc and Vallentin [TJHIS]). 

In fact, by using the stabilizer subgroup of the isometry group this method 
derives new positive semidefinite constraints which are stronger than constraints 
given by the Bochner - Schoenberg theorem (see PQ). In other words, it shows 
that constraints for {iy} and {ui} given by PD(M, {qi}) are stronger than 
constraints in PD(M, 0) for spherical, Hamming and some other spaces. 

Clearly, if Fi,F 2 <E PD(M, Q), then for any functions a\,a 2 in m variables 
we have 

F(t,u,v) = oti(u)a>i{v)Fi(t,u,v) + a 2 (u)a 2 (v)F 2 (t,u,v) € PD(M,Q). (1.1) 

In this paper we consider the class PD(§" _1 , Q). Namely, we show that this 
class can be generated by (1.1) using as a basis the polynomials G[™ , where 
m = card(Q) < n — 2. Actually, the class PD(S™~ 1 ,Q) becomes "more rich" 
whenever m is increasing. 

The paper is organized as follows: Section 2 introduces multivariate Gegen- 
bauer polynomials and its basic properties. Section 3 extends the Schoenberg 
theorem. Section 4 extends the Schoenberg theorem to Euclidean spaces. Sec- 
tion 5 considers the positive semidefinite constraints for the distance distribution 
given by Theorem 3.1. Section 6 discusses applications of the extended Bochner 
- Schoenberg theorem to spherical codes. 



2 Multivariate Gegenbauer polynomials 

In this section we introduce polynomials G^'™^ (t, u, v) in 2m + 1 variables with 
< m < n - 2, t e K, u = (m, . . . ,u m ), v = (vi, . . . ,v m ) G M m . 

2.1. Gegenbauer polynomials. There are many ways to define Gegenbauer 
(or ultraspherical) polynomials G^\t) (see [HI [SI EH E] ) • G^ are a special 
case of Jacobi polynomials Pjf with a — (3 — (n— 3)/2 and the normalization 
by (1) = 1. Also C?£ can be defined by the recurrence formula: 

ri-)-+ Mn) ( 2fc + " ~ 4) t ggi - (fc - 1) Gj-2 
G -1, G x G k k + n-3 ' 

Note that for any even k > (resp. odd) Gu (i) is even (resp. odd). 
Therefore, G 2 k and G 2 i+\ are orthogonal on [—1,1]. Moreover, all polynomials 
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G^ are orthogonal on [—1,1] with respect to the weight function (1 — t 2 )^ n 3 ^ 2 : 

1 G<™ } (t) G* n) (t) (1 - i 2 )("- 3 )/ 2 dt = 0, k £ £. (2.1) 
-l 

Recall the addition theorem for Gegenbauer polynomials: 
(cos 9\ cos 6 2 + sin #1 sin 9 2 cos </?) 

fe 

= $>„ fcs Gtt 2s) (cos0i)Gt+ 2s) (cos0 2 ) (sin^) 8 (sin0 2 ) s G^Hcos^), 

where c n ks are positive coefficients whose values of no concern here (see [5JIS])- 

2.2. Basic definitions. Let us denote by (u, v) the inner product of vectors 
u, v £ K"\ and by |v| we denote the norm of v (i.e. |v| 2 = (v,v)). 

Definition 2.1. Let < m < n - 2, < e M, u,v £ W n for m > 0, and 

u = v = for m — 0. Then the following polynomial in 2m + 1 variables of 
degree k in t is well defined: 

G { k n ' m) (t, u, v) := (1 - |u| 2 ) fc / 2 (1 - |v| 2 ) fc / 2 G { r m) ' f " ( "' v) I 



Vv(i-N 2 )(i-iv| 2 ); 

Note that by definition we have: G<"' 0) (i) = G^ 0) (t,0,0) = G^f). 
Let u = (u\, . . . , u n ) £ W 1 . For < m < n we denote by \i^> the vector 
... , u m ) in M m . Let u<°) = 0. For < m < n - 2, u, v £ E" put 

G<f' m) (t, u, v) := G^' m) (t, u(™) , v(™) ) . 

Definition 2.2. We define a matrix Z^(u, v) 0/ size ( m + d ) x ( m + d ) /or u £ 
R m , vel™ &y ^™(u,v) := (z7(u)) T z™(v), w/iere 

^rf — ("^ 1 ' ' ' ' > -^m) ■ — (1? ^1 ? • • • ) ) 5 X\X2 5 ■ • ■ ; %m—l ^rn 1 ^rn ' ' ' * ' ^m) 

is i/ie vector of monomials. 

Denote the inner product of matrices A — (<%) , B — (pij) both of size dx d 

by (A, B), i.e. (A, B) = Tr(AB) = £ n, ,b, r 

»>i=i 

Definition 2.3. Let /(u, v) &e a symmetric polynomial in u, v € R m , i.e. 
/(u, v) = /(v, u). We say £/ia£ /(u, v) is positive semidefinite and write f y if 
there are a symmetric matrix H >z and d > so £/ia£ /(u, v) = (if, Z™(u, v)) . 

Let polynomials f r >z /or aZZ r, / r — > / as r — > oo , where a function /(u, v) 
is symmetric in u, v and continuous on {|u| < 1, |v| < 1}. TTien we say £/iat 
/(u, v) is positive semidefinite and write f >z 0. 
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2.3. Orthogonality. Let isl, u = (ui,.. . ,u m ), v = (vi,. .. ,v m ) e 
Denote 

/ 1 ... Ml »i \ 



Q TO (t,u,v) := 



... 1 u m v m 
m ... u m 1 t 

\ Vl ... V m t 1 / 



Let 
Since 



D m := {(t,u,v) G R 2m+1 : Q m (t,u,v) h 0}. 



det(Q m (t,u,v)) = (1 - |u| 2 )(l - |v| 2 ) - (t- (u,v)) 2 , 
we have (t, u, v) g D m if and only if 

(t-(u,v» 2 <(l-|u| 2 )(l-|v| 2 ), |u|<l. 

Denote (det(Q TO (i,u,v)))("-™- 3 )/ 2 by p„ ;TO (t,u,v), i.e. 

p„ ;ro (i,u,v) := ((1 - |u| 2 )(l - |v| 2 ) - (t- (u,v)) 2 ) (n - m - 3)/2 . 

Theorem 2.1. Let < m < n — 2. Let q(u,v) be any continuous function on 
{(u, v) g R 2m : |u| < 1, |v| < 1}. If k^e, then 

/ G^ l,m) (t, u, v) G { P' m) (t, u, v) q(u, v) p„ !m (t, u, v) dt du dv = 0. 
JD m 

Proof. Let 

t = S ((l-|u| 2 )(l-|v| 2 )) 1 /2 + (u , v) . 

In the variables s, u, v we have 

D m = {(s,u,v) : -1 < a < 1, |u| < 1, |v| < 1}. 
The Jacobian of this change of variables is ((1 — |u| 2 )(l — | v| 2 ) ) . Then 

1=1 G^ m) (t,u,v)G ( t m \t,u, v) g(u, v) p n , n (t, u, v) dt du dv 

= I J' Gi n ~ m) ( S ) G ( p- m) ( S ) (1 - s 2 ) (n-m-3)/2 ^ 

Io = I I ((l-|u| 2 )(l-|v| 2 )) ( "- ro+fe - 2)/ %(u,v)durfv. 

J |u|<l J\v\<l 

Thus, (2.1) yields 7 = 0. □ 

Let ei, . . . , e„ be an orthonormal basis of R™, and let (x\, . . . , x n ) be the 
coordinates expression of a point x G R™ in this basis. Then x {m ^ as well as 
G ( "' m \(x,y),x,y) = G { ™' m \{x, y), x^ , y^) for x,y G R" are well defined. 
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Lemma 2.1. For any continuous F on D m and < m < n — 2 we have 
f F((x,v),xW,yW)du, n (x)dLJ n (y) 

J(S«-1)2 

/ F(t,u,v)r m (u) r m (v) p n , m (t, u, v) rfi du dv, 

./-Dm 

w/iere r m (u) = n£Li (-*■ — |u^) | 2 ) !/ 2 f or m > 1, ro = r\ = 1, and to n is the 
surface area o/S™ -1 for the standard measure ckv n . 

Proof. Let x,y £ S"" 1 . Let t = (x,y), u = x<- m \ v = y( m ). Let 

a = — — , b = -r-—-, s = (a, 6), where h(u) = y/l - |u| 2 . 
/i(u) a(v) 

Then 

y), x^ , yW) = F(ft(u) h(v) (a, b) + (u, v), u, v) = F(t, u, v). 

Therefore, we can change the variables (x,y) — > (a, 6, u,v) — > (s,u, v) — > 
(i,u,v). 

We obviously have |u| < 1, |v| < 1, -1 < s < 1. Clearly, \a\ = \b\ = 1. 
Then (a, 6) e (S"-™- 1 ) 2 . Since |f - (u, v) | < h(u) h(v) we have (t, u, v) £ D m . 
Consider (x, y) — > (a, 6, u, v). The Jacobian is defined by 

doj n (x) = fc„, m (u) dudu;„_ m (a), 

m— 1 

fen, m (u) = (1 - |u| 2 )(" — 2 )/ 2 [] (1 lu^l 2 ) 1 / 2 . 

i=0 

Combining the volume of {(a, 6) £ (§" _m_1 ) 2 : (a,b) = s} and the Jacobian of 
(a, b, u, v) — > (s, u, v) we get: w n _ TO a>„_ m _i p„ >m (s, 0, 0). The Jacobian of the 
changing (s,u, v) — > (£, u, v) is (h(u) /i(v)) -1 . That completes the proof. □ 

Combining Lemma 2.1 and Theorem 2.1 we get: 

Corollary 2.1. Let < m < n — 2. Lei /(u,v) &e any continuous function on 
{(u, v) G R 2m : |u| < 1, |v| < 1}. If k^£, then 

f G^ m \(x, y),x, y) G^ m \(x, y),x, y) f(x™ , y^) dco n (x) dcu n (y) = 0. 
2.4. The addition theorem. 

Theorem 2.2 (The addition theorem). Let 1 < m < n — 2. Then 
G^-'Ht, u, v) = £ C fc ":"» CH(v) G<"-) (i, u, v), 
w/iere C^' m (u) = C^' m (ui, . . . , u m ) is a polynomial in u\, . . . , u m of degree d. 
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Proof. Suppose t, u = (tii, . . . , it m ), v = (v\, . . . , w m ) are such that 
|u|<l, |v|<l, (t-(u,v» 2 <(l-|u| 2 )(l-|v| 2 ). 
Then 6\,&2,<p € (0, 7r) are uniquely defined by the following equations: 



cost/i = — ; =, cos #2 = 

l-uf-.-.-w^i v - 1 - r'f - <„,. 

t-(u,v) 



COS(,5 



V(i-K)(i-| v p) 

The addition theorem for Gegenbauer polynomials yields: 

Gt m - 1] (t, u, v) = £ C£_*» C£_*» G<"'™) (i, u, v), (2.2) 

s=0 



s-m,m 1 \ 1 k — s ^i(n+2s) / U m \ H. 2 

C fe i s (u) := V^Iw G^_ s J, W = ^l-u\ 



*ra-l- 



It's easy to see that C^i™(iti, . . . , u m ) is a polynomial of degree k — s. Thus, 
(2.2) holds for all t,u,v. □ 

Let b s (u, v) = C fc "l™(u) C*fcl™(v). Clearly, b s £ 0. Then Theorem 2.2 yields 

Corollary 2.2. Let < m < £ < n - 2. Let u, v e M™. T/ien 

Gl"' m) (t , u, v) = ^ /. ( U M , v M ) G («.0 (t , u , v) 

f s £ for all < s < k. 

3 An extension of the Schoenberg theorem 

3.1. Schoenberg's theorem. Let p±,...,p r be points in §™ _1 , and let 
aij . . . , a r be any real numbers. Then 

< \aipi + . . . + a r p r \ = (pi,pj)a,iaj, 

id 

or equivalently the Gram matrix (jflijPj)^ is positive semidefinite. 

Schoenberg [17] extended this property to Gegenbauer polynomials G£ . He 
proved that for any finite X — . . . ,p r } C § n_1 the matrix ^G^\{Pi, Pj))^ 
is positive semidefinite. 

Schoenberg proved also that the converse holds: if f(t) is a real polynomial 
and for any finite X C the matrix (f({pi,Pj))) £ 0, then f(t) is a linear 

combination of G^ 1 (t) with nonnegative coefficients. 

3.2. An extension of the direct Schoenberg theorem. 
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Theorem 3.1. Let e\, . . . , e n be an orthonormal basis o/R™, and let pi, . . . ,p r 

be points in § n . Then for any k > and < m < n — 2 the matrix 

■ /'.; •/'-•/'. )) is positive semidefinite. 

Proof. Actually, this theorem is a simple consequence of the Schoenberg theo- 
rem. Indeed, let Vj — (pi, e\) e\ + . . . + (p^, e m ) e m , and let x,; = pi — v^. Then 
Xi is a vector in R ra_m with the basis e m +i, . . . , e„. Note that | | 2 + |xi| 2 = 1. 
Let y; = Xi/|xi| for |xj| > 0. In the case |x^| = put = e„. Then 

Recall the Schur theorem: If A ~ Wj) ^ 0, B = (by) ^ 0, t/ien C = 
(a^&y) h 0. Let ftj = | Xi | fc = (1 - |v,| 2 ) fe / 2 , h = {h ll ...,h r ), A = h T h. 

Clearly, A>Q. The Schoenberg theorem yields: B = (Gjf~ m) ((y i; y^-))) ^ 0. 

It is easy to see that G£' (£y,Vj,Vj) = ay&ij, where ty := (pi,Pj)- Thus, 

C=(4^(^, Vi , Vi ))h0. ^ 

Using the addition theorem it is not hard to give a direct proof. Note that 

(2) (2) 

G k is the Chebyshev polynomial of the first kind, i.e. G k (cos (9) = cos(kd). It 

is easy to see that for any <pi, . . . , tp r the matrix ^cos(< ( 5 i — (pj)j b 0. From this 

follows that Schoenberg's theorem holds for n = 2. Therefore, we have proved 
the theorem for m = n — 2. Put I = n — 2 in Corollary 2.2. That yields the 
theorem for all < m < n — 2. □ 

Remark 3.1. We see that the second proof (as well as Schoenberg's original 
proof in [17] ) is based on the addition theorem for Gegenbauer polynomials. 
There exists another proof of Schoenberg's theorem which is using the addition 
theorem for spherical harmonics (see, for instance, |16j). It is possible (see [I] 
for the case m — 1) to derive Theorem 3.1 from this theorem. However, this 
proof looks more complicated than our proof. 

Remark 3.2. Bachoc and Vallentin [T] derived new upper bounds for spherical 
codes based on positive semidefinite constraints that are given in [TJ Corollary 
3.4]. In fact, this corollary easily follows from Theorem 3.1 with m = 1. 

Indeed, let us fix some £, 1 < I < r, and take an orthonormal basis ei, . . . , e n 
of M. n with ei = pe- Then for m=lwe have = tn. Let 

( A i)ij = G^ A) (fij,tu,tjt), 0<i< r ,l<j<r, 

Y e = WA e W T , {W) l3 = A 4 G|" +2fc) (^), < * < d, 1 < j < r. 

Theorem 3.1 yields: A t > 0. Then Y t > and Yi + . . . + Y r h 0. These 
constraints are equivalent to the constraints in [TJ Corollary 3.4]. 

Corollary 3.1. Let < m < n — 2, d > 0. Let e\, . . . , e n be an orthonormal 
basis o/R™. Let a polynomial F(t,u,v) can be represented in the form 

d 

F(;,u,v)=^/ fc (u,v)G^ l) (i,u,v), 

k=0 
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where fk >z for all < k < d. Then for any points p\ , . . . , p r in S" 1 
the matrix ^F((pj,pj),^ m \pj m ^ is positive semidefinite. In other words, 
F ePB(S n -\{ ei ,...,e m }). 

Proof. Let /(u,v) h 0. Then /(u,v) = (H,Z^(u,v)) with H > 0. From an 
eigenvalue factorization of if it follows that there exist polynomials hi (u) such 

that /(u,v) - £ifci(u)Mv). Therefore, (/(p| m) ,p< m) )) b 0. 
We have A fe = (/k(pj m) ,pj m) )) >: 0. Theorem 3.1 yields 

B k =(G^((p^ P] ),pt\p { r ) ))hO. 
Let (Ck)ij = (Ak)ij (P>k)ij- The Schur theorem implies: Ck h 0. Thus, 

v 7 fc=0 

3.3. An extension of the converse Schoenberg theorem. 

Theorem 3.2. Let < m < n — 2. Let e\, . . . , e n G M n &e an orthonormal basis. 
Let F(t, u, v) fee a polynomial in t and a symmetric polynomial in u, v € R m . 
Suppose F e PD(S"-\ {ei, . . . , e m }). T/ien 

d 

F(t, u, v) = ^ / fc (u, v) G<"< m) (i, u, v) (3.1) 
with f k h0 for all < k < d = deg t (F). 

Proof. The polynomial G k n ' m \t, u, v) has degree fc in the variable t, so that 
F has a unique expression in the form (3.1), where /fe(u, v) is a symmetric 
polynomial in u, v of degree du in u. Then /fc(u, v) = (Fk, Z™ (u, v)}, where 
Fk is a symmetric matrix. Now we prove that Fk >z 0. 

Let f(t, u, v), g(i, u, v) be continuous functions on _D m . Denote 

\f,g\ := I f((x,y),x^\y^)g((x,y),x^\y^)diJ n (x)du; n (y). 

For any continuous function h(u, v) Corollary 2.1 yields: 

{FhG k n ' m) } = {f k G k n ' m \hG k n ' m) } = [f k K (G k n ' m) ) 2 } • (3.2) 

Let /(u,v), g(u, v) be continuous functions on {|u| < 1, |v| < 1}. Denote 

[f,g] :={fg,(G k n ^) 2 }. 

Clearly, [•, •] is an inner product. Let a\, a.2, ... be an orthonormal basis for [•, •] 
in the space of real polynomials R[u, v]. We observe that Fk defines a quadratic 
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form on R[u,v]. Let us denote by F k the matrix expression of this quadratic 
form in the basis {at}. Then there exists a matrix A such that F k — A T F k A. 
Therefore, F k h yields F k h 0. 

Let h(u,v) = (H,Z^ k (u,v)). It's easy to see that [f k ,h] = (F k , H) . 

Note that the sum of all entries of a positive semidefinite matrix is nonnega- 
tive. This implies: if / G PD(S n -\ E m ), E m := {ex, . . . ,e m }, then {/, l} > 0. 
Moreover, the Schur theorem yields: if / and g £ PD(§™ -1 , E m ), then fg G 
PD(S n - 1 ,E m ). Therefore, {f,g} > 0. Thus, for any f,g G PD(§ n -\£: m ) we 
have [f,g] > 0. 

Let H y 0. Then ft(u,v) G PD(§™~ 1 , E m ). We have F G PD(S n -\ E m ). 
Then (3.2) yields [f k ,h] = (F k , H) > for any H h 0. Thus, F k b and 
F fc b 0. □ 

Remark 3.3. In 2, Proposition 4.12] a statement that is equivalent to Theorem 
3.2 was proven with m = 1. However, in [21 Proposition 4.12] the formula (3.1) 
is written in terms of an orthogonal basis for {•,•}• 

3.4. The class PD(§"- 1 ,Q). If f r G PD(M, Q), f r -> / as r oo, and 
f(t,u, v) is continuous, then also / G PD(M, Q). Therefore, Corollary 3.1 and 
Theorem 3.2 imply 

Theorem 3.3. Let < m < n — 2. Let ei, . . . , e„ be an orthonormal basis of 
IT . Then F £ PD(§™~ 1 , {e x , . . . , e m }) i/ and only if 

CO 

L(t,u,v)=^/ fc (u,v)G^ m) (t,u,v), 

fc=0 

where for any k > a function f k (u,v) is positive semidefinite. 

Now we consider PD(§"~ 1 , Q) for any Q = {q 1: q m } C S™" 1 with < 
m < n — 2. Theorem 3.3 yields 

Corollary 3.2. Let < m < n - 2. Let Q = {qi,...,q m } C S" _1 wift 
rank(Q) = m. Let e\, . . . , e m be an orthonormal basis of the linear space with 
the basis q%, . . . , q m , and let Lq denotes the linear transformation of coordinates. 
Then F G PD(S"" 1 , Q) if and only if 

OC 

F(t, u, v) = /*(u, v) G k n ' m) (t, L Q (u), L Q (v)), 

fe=0 

where /fc(u, v)^0 /or all k > 0. 

Remark 3.4. It is not hard to describe PD(§ ,1_1 , Q) with rank(Q) > n — 1. 

First, consider the case Q = {ei,...,e„_i}, where {e{\ is an orthonormal 
basis. Let p x , . . . ,p r G S n_1 , U; = (pi, a) e x + . . . + (pi, e n -i) e„_i, and Uj = 
(pi,Pj). Obviously, Uj - (Uj,Uj) = WiW 3 , Wi = (pi,e n ). Then 

^-(u,,^)! =w T w^0, w = (wi,...,w r ). (3.3) 
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Since w\ = 1 — | | 2 , we have for alH = 1, . . . , r, j — 1, . . . ,r: 

(^-(u l; u,)) 2 HlHud 2 )(l-N 2 ). (3.4) 

Let fla(<,u,v) := (i - (u,v)) 2 - (1 - |u| 2 )(l - |v| 2 ). Since H 2 (t ij ,n h u j ) = 0, 
we have F G PD(S n_1 , {ei, . . . , e„_i}) if and only if 

F(t, u, v) = f (u, v) + /i (u, v) (t - (u, v» + ii(t, u, v) H 2 (t, u, v), 

where /o _: 0, /i 0, and R(t, u, v) is any continuous function on D m . 

Let Q — {qi, . . . , g n -i} C with rank(Q) = n — 1. Let 1/ be a linear 

transformation of the basis qi, . . . , q n -i to an orthonormal basis e±, . . . , e n -\ of 
R™- 1 . Then F £ PD(§ n_1 ,Q) if and only if 

F(t, u, v) = / (u, v) + /x(u, v) (t - (L(u), L(v)» + R(t, u, v) H 2 (t, L(u), L(v)), 

with f h0, fit 0, and any R(t, L(u), L(v)) G C(D ro ). 

In the case rank(Q = {91, . . . , q m }) > n — 1 consider all Qi C Q with 
rank(Qi) = n — 1. Let a linear transformation Li : M m — + R" _1 is defined by 
Li{q) = q if q G Qi, = if q G Q and rank(Qi U {g}) = n. Denote by N 

the number of distinct Qi. It is not hard to prove that F G PD(S" _1 , Q) iff 

N 

F(t, u, v) = £ /i(u, v) f}(t, L<(u), Li(v)), 

i=l 

where for all t: /i(u,v) >r 0, Fi(i, L;(u), Lj(v)) G PD(S"- 1 ,Q i ). 

4 Positive definite functions in W l 

Direct extensions of the Bochner - Schoenberg theorem and finding bounds 
on sphere packings in R™ are not so straightforward because this space is not 
compact. Different indirect ways of deriving bounds on sphere packings in R n 
were suggested in the literature [El EJ U\ ■ 

Here we note that the multivariate Gegenbauer polynomials defined above 
enable one to define p.d. functions in R™ as follows: 

Hl n > m \t,x,y,u,v) := (xy) k / 2 (f,u', v'), 

where < m < n — 2, t, x, y G R, u, v G R m for m > and u = v = for 
m = 0, and 

t' = t/y/xy, u' = u/y/x, v' = v/v/y. 

The positive semidefiniteness of the polynomials Q^ ,m ^ (Theorem 3.1) implies 
the following result. 

Theorem 4.1. Let e%, . . . , e n G R™ be an orthonormal basis, and let p\, . . . 

be points in R n . Then for any k > and < m < n — 2 the matrix (jjij) , where 

9ij = H i ; (<Pi,Pj), \Pt\ 2 , \Pj\ 2 ,Pi,Pj), 

is positive semidefinite. 
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Remark 4.1. It is not hard to find Euclidean analogs of Theorems 3.2, 3.3 and 
Corollary 3.2 

Theorem 4.1 gives a family of positive-semidefmite constraints for distance 
distributions of points in Euclidean spaces. For instance, consider the simple 
case of m = 0. Now for any matrix A = (ay) of size N x N we have a matrix 

h£\A) which is defined by 

Corollary 4.1. If A is a symmetric positive semidefinite matrix with rank(A) < 
n, then for any positive integer k we have 

Hi n \A)yo. 

Proof. It is a well known fact: a symmetric matrix A of size N x N is the Gram 
matrix of N vectors in R n if and only if A >z and vank(A) < n. Therefore, 
there are vectors p\, . . . ,pn in R n such that A = ({Pi,Pj})- Then 



(Hi n \A)) i .=Hi nfi \(p i ,p j ),\p i \\\p j \^0,0). 

r(n) 



Thus Theorem 4.1 yields Hf> (A) h 0. □ 

For small k it is not hard to give explicit expressions for H^ 1 (A) . Clearly, 
H[ n) (A) = A. Since G<?\t) = (nt 2 - l)/(n- 1), we have 

V 2 V )) tJ n l 

Then 



flj-) (j l)-!^.-« T » 



n-l ' 
where 

(A 2 )ij := a?-, a := (an,a 2 2, • ■ • , aArAr). 
So if A y and rank(yl) < n, then 

nj4 2 >: a T a. 



5 Positive semidefinite constraints 

In this section we consider positive-semidefinite constraints that are given by 
Theorem 3.1. It is a natural to ask which constraints are stronger than the 
others? We will show that if mi > m 2 , then the constraints for m = m\ imply 
the constraints for m = m 2 . 

Let T = (tij) be a symmetric matrix of size rxr with — 1 < tij < 1, ta = 1, 
and let U — {u%j) be a matrix of size r x (n— 1). Let < 1 for alH = 1, . . . , r, 
where Ui := (un, . . . , Ui, n -i)- Then we say that a pair (T, U)" is feasible. 



11 



Proposition 5.1. Let < to < £ < n - 2, d > 0. Let {T,U) r r l be a feasible 
pair. Suppose 

(Gj^ty.Ui.u,-)) hO for k = l,...,d. 

Then 

Ol' m {T, U) := (Gjj. n,m) (t« , Ui , Uj)) ^0 /or a// fc = 1, . . . ,d. 

Proof. Corollary 2.2 (see also the proof of Corollary 3.1) yields: G"' m (T, U) is 
a sum of positive semidefinite matrices. □ 

In Remark 3.4 we considered the case to = n — 1. Now we show that the 
constraints (3.3), (3.4) are strong enough. 

Proposition 5.2. Let (T, [/)" &e a feasible pair. Suppose 

T h ((uj,Uj}), 
cmd /or all i = 1, . . . ,r, j = 1, . . . , r, we /lave 

(^•-(u^^a-KI'Xi-M 2 ). 

ITien i/iere are points p\, . . . ,p r , e\, . . . , e n -i in §™ _1 with (ei,ej) — Sij such 
that tij = (pi,Pj), u lk = (pi,e k ) for all i,j,k. 

Proof. Consider a symmetric matrix Xq — (xij) of size £ x £, £ = r + n — 1, 
that is defined by: = t^ for 1 < i < r, 1 < j < r; = 8ij for r < i < 
£, r < j < £; x^ = Ui s , s = j — r + n — 1, for 1 < i < r, r < j < £; and 
x^ = Xji for r < i < £, 1 < j < r. 

It is easy to sec that (3.4) yields rank(Xo) < n. On the other hand, it follows 
from (3.3) that X n y 0. Therefore, there are vectors qi, . . . , q r + n -i in M™ such 
that X n = ((qi,qj))- Since xu = 1, we have \qi\ — 1, i.e. qi e § rl_1 . Denote 
Pi = qi for i = 1, . . . , r, and = q i+r for i = 1, . . . , n — 1. □ 

In particular, we obtain that (3.3), (3.4) imply G^' m (T, {/) y for all to : 

< to < n- 2. 

Actually, G^' m (T,U) y gives constraints only for T and u^ m \ For in- 
stance, if to = 0, then we just have | | < 1. Now using Theorem 3.1 we 
improve constraints for U. 

Let ei, . . . , e n be an orthonormal basis of K™, and let pi, . . . ,p r be points 
in S n_1 . Denote t tj = {p t ,p 3 ), 1 < i < r, 1 < j < r, u ik = (pi,e k ), where 

1 < k < n — 1 . Then we have the matrices T = (t^ ) and U = (iiij ) . 

Consider points {qi, . . . ,q e } = {pi, . . . ,p r ,e m+ i, . . . ,e„_i} in S n_1 , where 
< to < n — 2, £ = r + n — to — 1. Then ccy = (qi,qj) and = (qi,e k ), 1 < 
k < m, define matrices 

X m = X m (T, (7) := (a:«) , V m = V m (T, U) := (« y ) . 

Here y m is well defined for to > 0. Put Vq = 0. 
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Let v, := (vn, . . . , i>, m ), where 1 < i < I. If we apply Theorem 3.1 for the 
points {qi, . . . , qi}, then for any k > we get 

G^ ro (X m ,y m ) := (G^^-.Vi.v,-)) >: 0. 

Let (T, [/)" be a feasible pair. Then X TO (T, J7) and V m (T, U) are well defined. 
It is clear that 

G n k ' m {X m {T, U), V m (T, U)) h yields G n k ' m (T, U) h 0. 

Denote by A^'™, < m < n — 2, d > 0, r>0, the set of all feasible pairs 
(T, (7)™ such that' 

Gl m (X m (T,U),V m (T,U))hO for fc=l,...,d. 

Proposition 5.3. A^ l ~ 2 c . . . C A^ 1 c A£°. 

Proof. It is well known fact: a symmetric matrix A = (a^) o/ size I x f wii/i 
a« > is positive semidefinite if and only if B = (pij) y 0, where 

6 .. = a .._£i£M ) l < i,i < 1. (5.1) 

Let (T, <7) e AJ™ +1 , where < m < n - 3. Denote 

A = Gr i (A%„(T,;y),v m (T,[/)). 

Then from the addition theorem (Theorem 2.2) we have: 

<>, >>.,,• fc«.*J =CfciT H1 (vi)C^lT H1 (v J -)G^+ 1 )( a: «,v i ,v J -). 

Note that c n ko = 1, where c n fc s are coefficients in the addition theorem for 
Gegenbauer polynomials (see subsection 2.1). Indeed, put 9\ = 62 = <p = 0. 
Then we have 1 = G^ l) (l) = c nk0 (G k n) (l)) 2 G^™ _1 ^(1) - c„ fc0 . This yields: 

ho,i,j = -, 1<«,J<^-1- 

au 

Let B is defined by (5.1). Then bij = Y^=i h a ,i,j- Therefore, B — J2 S=1 H s , 
where H s — (h s ,i,j)- Since H s y 0, we have B y 0, ^4 y 0. Thus we have 

j a n.m+1 ~- * n,m i — i 

proved: A d ' r C A d ; r . □ 

We denote by A™ the set of all feasible pairs (T, t/)™ that satisfy the as- 
sumptions of Proposition 5.2 (i.e. A" consists of pairs (T, £/)" that satisfy 
(3.3), (3.4)). 

Note that for 1 < m < n— 2 we have (" ) choices for u^" 1 ). In other words, 
not only ei, . . . , e m , but any m vectors from {ei, . . . , e„_i} can be chosen as a 
basis of W n . Denote by SAJJ'™ the intersection of all corresponding spaces AJ™. 
It is not hard to see that Propositions 5.2, 5.3 imply 

Proposition 5.4. A? c SA n d '^ 2 C...C SA# c A™' r °. 
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6 Upper bounds for spherical codes 



In this section we set up upper bounds for spherical codes which are based on 
multivariate p.d. functions. These bounds extend the famous Delsarte's bound. 
Note that for the case m = 1 this bound is the Bachoc - Vallentin bound [3]. 

Definition 6.1. Consider a vector J = (Jx, . . . ,jd). Split the set of numbers 
{ji, . . . , jd\ into maximal subsets Ii, . . . , Ik with equal elements. That means, 
if I r = {j Tl , ■■■,jr,}i then j ri = . . . = j rs = a r and all other ji ^ a r . Without 
loss of generality it can be assumed that i\ = > . . . > ip. = \Ik\ > 0. (Note 
that we have %\ + . . . + i& = d.) Denote by ip(J) the vector u> — (ix, . . . ,ife). 
Let 

W d ■■= {u = (h, ...,ik)'-ii + --- + ik = d, ii > . . . > i k > 0, h, .. . ,i k € Z}. 
Let lu £ Wd- Denote 

UN) := #{ J = (,7i, ■ • • e {1, • • • , N} d : = w}, 

UN) 



N 



It is not hard to see that g w ( N ) is a polynomial of degree d — 1 for to € Wd 
and 

qUN) = N d -\ 

u)£W d 

Definition 6.2. For any vector x = {xij} with 1 < i < j < d denote by A(x) 
a symmetric d x d matrix (dfj) with all an = 1 and aji = aij = x^, i < j. 
Let < 9 < 7r and 

X(9) := {x = {x^} : x^ G [— l,cos6*] or x^ — I, 1 < i < j < d}. 

Now for any x = {x^} G X{9) we define a vector J(x) = . . . ,jd) such 
that jk — k if there are no i < k with Xik = 1, otherwise jk — i, where i is the 
minimum index with Xik = 1 ■ 

Let uj G Wd- Denote 

D u (6) := {x G X{9) : i>(J(x)) = uj and A(x) b 0}. 
Let /(x) be a real function in x, and let 

B U (6J):= sup /(x). 

xED„(fl) 

Note that the assumption A(x) y implies existence of unit vectors pi, . . . ,pd 
such that ^4(x) is the Gram matrix of these vectors, i.e. Xij — (pi,pj). More- 
over, if Xij = 1, then pi = pj. In particular, D^d){9) — {(1, . . . , 1)} and therefore 
B {d) {6J) = /(!,..., 1). 
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Definition 6.3. Let x = {:%}, where 1 < i < j < rn + 2 < n, and let A(x) >z 0. 
Then there exist P = {p\, . . . ,p m +%} C § n_1 such that Xij — (pi,pj). Let 
F(x) be a continuous function in x with Ffeki) = -F'(x) for all x/^ that can be 
obtained by interchanging two points pk and pg in P. We say that F(x) £ PD^ 
if for all x with A(x) >z we have F(xi2, Ui, U2) G PD(§ ,1_1 , Q(x)), where 
u ; = (x i3 , • • ■ , x iiTn+2 ), Q(x) = {p 3 , . . . ,p m+2 }, and F(xi 2 ,u 1 ,u 2 ) = F(x). 

For the classical case m = Schoenberg's theorem says that / € PDq ^ an d 

only if f(t) = J2k /fe G fe n) (*) with a11 fk > 0. Theorem 3.3 (see also ED |3]) 
yields F(x) e PD™ if and only if 

F(xi2, X 13 , X 23 ) = ^ fk( X 12,Xi 3 ,X 2 3) G < £ l,1 \xi2, X 13 , X 23 ), 
k 

where ft > for all k, and 

F(xi2, X 13 , X 23 ) = F(x 13 , Xi2,X 2 3) = F(x 23 ,xi 3 , X12). 

Using Corollary 3.2 it is possible to describe the class of functions PDJ^, m > 2. 

Let C be an iV-element subset of the unit sphere § n_1 C 1". It is called 
an (n, N, 9) spherical code if every pair of distinct points (c, c') of C have inner 
product (c, c') at most cos 6*. 

Theorem 6.1. Let f > 0, < m < n - 2, and F(x) = /(x) - / € PD^. 
Then an (n, N, 9) spherical code satisfies 

foN m+1 < £ B u (0,f)qUN). 

Proof. Let C be an (n, N, 9) spherical code. Define 

S = ]T f({(ci, Cj )}), d = m + 2. 
c=(c u ...,c d )ec d 

Then 

S = E E /({< c i» c i)»< E BU0J)UN). 
On the other hand, since F € PD^ we have 

J2 F({(c i ,c,-)})>0. 

Thus 

5= E (/o + n{(^c,)}))>/ iV d . 
cec d 

□ 
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It is easy to see for m = that q( 2 )(N) = 1, q(i t i)(N) = N — 1, and 
B(2)(0,f) — /(l)- Therefore, from Theorem 6.1 we have 

/oiV</(l) + B(ia)(0,/)(^-l). 

Suppose B(i,i)(0,/) < 0, i.e. /(t) < for all f € [-l,cos0]. Thus for (n,N,9) 
spherical code we obtain 

~ fo 

This upper bound is called Delsarte's bound. 

The Bachoc-Vallentin bound [H Theorem 4.1] is the bound in Theorem 6.1 
for to = 1 and B {1 1 1} (0, f) < 0. Indeed, let B (2 1} (0, /) < B. Since <j (3) (N) = 
1, 9(2,1) (AO = 3(V- 1), and B (3) (0,/) = /(l, 1,1), we have 

foN 2 < f (1,1,1) + 3(N-1)B. 

Let us consider Theorem 6.1 also for the case m = 2 with i?n 1,1,1) (6, f) < 0. 
Let 5(3,1)^,/) < Bi, B (2>2) (6,f) < B 2 , and B (2tlA) (6,f) < B 3 . Then 

/oiV 3 < /(1, 1, 1, 1, 1, 1) + 4(iV - l)Bi + 3(iV - 1)B 2 + 6(N - 1)(N - 2)B 3 . 

Let /(x) be a polynomial of degree d. Then the assumptions in Theorem 
6.1 can be written as positive semidefinite constraints for the coefficients of F 
(see for details [TJ [31 [101 dH] ) • Actually, the bound given by Theorem 6.1 can 
be obtained as a solution of an SDP (semidefinite programming) optimization 
problem. In [TJ [2] using numerical solutions of the SDP problem for the case 
to = 1 has obtained new upper bounds for the kissing numbers and for the 
one-sided kissing numbers in several dimensions n < 10. 

However, the dimension of the corresponding SDP problem is growth so fast 
whenever d and to are increasing that this problem can be treated numerically 
only for relatively small d and small to. It is an interesting problem to find 
(explicitly) suitable polynomials F for Theorem 6.1 and using it to obtain new 
bounds for spherical codes. 
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